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Abstract. In this paper, a discrete state transition algorithm is intro- 
duced to solve a multiobjective single machine job shop scheduling prob- 
lem. In the proposed approach, a non-dominated sort technique is used 
to select the best from a candidate state set, and a Pareto archived strat- 
egy is adopted to keep all the non-dominated solutions. Compared with 
the enumeration and other heuristics, experimental results have demon- 
strated the effectiveness of the multiobjective state transition algorithm. 
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1 Introduction 

The multiobjective optimization is encountered in many real world applications 
PP. For a specific policy, the decision maker may find it advantageous for one 
goal but disadvantageous for others. A traditional way to deal with this issue is 
to impose a priori preference reflecting the relative importance of different ob- 
jectives; however, the final solution just indicates a decision maker's satisfaction, 
and it might be dissatisfactory for other decision makers. 

To ameliorate the problem, the concept of Pareto optimality and other rele- 
vant concepts are introduced. These are defined as follows: 

1) Pareto dominance: A feasible solution x = (x\, ■ ■ ■ ,x n ) is said to Pareto 
dominate another feasible solution y = • • • , y n ), denoted as x -< y, if 

/i(aO</<(tf),V»e{l,--. ,k}, and 3j E {1, • • • , k}, f 3 (x) < f 3 (y), (1) 

where, fi(x) is the ith objective function, k is the number of objectives. 

2) Pareto optimality: A feasible solution x* is said to be Pareto optimal if 
and only if 

-dx G S,x -< x*, (2) 



* Corresponding author for this paper. This research work is conducted between 
Deakin University and Ballarat University under the Collaboration Research Net- 
work (CRN) initiative. The problem studied in this paper is related to the Australian 
Research Council (ARC) linkage project number LP0991175. 
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where, S is the feasible space. 

3) Pareto optimal set: The Pareto optimal set, denoted as P* , is defined by 

P* = {x* e S\-Bx eS,x-< x*}. (3) 

4) Pareto front: The Pareto front, denoted as Pf* , is defined by 

Pf* = {(fi{x*),--- ,f k (x*))\x* £P*}. (4) 

The introduction of Pareto optimality allows us to find a set of Pareto optimal 
solutions simultaneously, independent of the decision maker's priori preference. 

In the past few decades, evolutionary-based and nature-inspired multiobjec- 
tive optimization techniques have drawn considerable attention for scheduling 
problems 2 3 4 5 6 7 . In this paper, we introduce a recently new heuristics called 
state transition algorithm [8 9110111] as the basic search engine for the multiob- 
jective optimization. A non-dominated sort approach is used to select the best 
from a candidate state set, and the best state is stored using a Pareto archive 
strategy. Experimental results have testified the effectiveness of the proposed 
algorithm. 

2 Problem Description 

In the field of joinery manufacturing, jobs with similar materials can be scheduled 
together to minimize the amount of materials used; therefore, reducing the cost. 

For example, based on the cost savings matrix shown in Table [TJ pairing Jobl 
and Job2 will provide saving in the cost equivalent to 4 units. 



Table 1. The cost savings matrix for 5 jobs having the same material 





Jobl 


Job2 


Job3 


Job4 


Job5 


Jobl 





4 


2.64 


4.08 


3.9 


Job2 


4 





3.64 


4.72 


4.23 


Job3 


2.64 


3.64 





2.65 


2.87 


Job4 


4.08 


4.72 


2.65 





3.84 


Job5 


3.9 


4.23 


2.87 


3.84 






Additionally, based on the jobs' processing times and due dates as shown in 
Table [21 and for any given sequence and pairing of jobs, not only the total cost 
saving C is affected but also the total tardiness time T, which is calculated as: 

n 

T = ^2max{0,c j -d j } (5) 
where, Cj and dj are the completion time and the due time of job j, respectively. 
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Table 2. Due dates and processing times for a set of 5 jobs 



Job 


Due Date(days) 


Processing Time (hours) 


Jobl 


8 


17:40 


Job2 


2 


24:00 


Job3 


11 


19:20 


Job4 


3 


25:00 


Job5 


3 


14:40 



a Number of operational hours = 8 hours per day 



The goal of this paper is to determine the optimal sequence with pairing, in 
order to maximize the total cost savings and minimize the total tardiness time. 

It is obvious that finding the permutation of the sequence {1, 2, • • • , n} with 
pairing becomes a solution to the multiobjective single machine scheduling prob- 
lem; however, not without the necessity to discuss the number of pairs for any 
fixed sequence of jobs. 

Given a sequence s = (1,2, •• • ,n), for n = 3, we have 2 possible pairing 
options (l-2)-3 and l-(2-3); for n = 4, we have 2 possible pairing options (1-2)- 
(3-4) and l-(2-3)-4, as pairing options (l-2)-3-4 and l-2-(3-4) are discarded; for 
n = 5, we have 3 possible options (l-2)-(3-4)-5, (l-2)-3-(4-5) and l-(2-3)-(4-5), 
as options (l-2)-3-4-5, l-2-(3-4)-5, l-2-3-(4-5) and l-(2-3)-4-5 are discarded. 

If Pl(n) denotes the number of pairs with the first two jobs pairing, and 
P2(n) denotes the complement of Pl(n), then we have the following theorem: 

Theorem 1 

Pl(n + 1) = P(n - 1), P2(n + 1) = Pl(n),n > 3 (6) 

where, P(n) = Pl(n) + P2(n) is the total number of pairs. For example, -PI (2) = 
1,P2(2) = 0, PI (3) = 1,P2(3) = 1, Pl(4) = 1,P2(4) = 1, Pl(5) = 2,P2(5) = 
1, we have Pl(4) = P(2),P2(4) = Pl(3), PI (5) = P(3),P2(5) = PI (4). 

Fig. [TJ shows the growth trend of the number of pairs with the sequence 
size; however, only small size job scheduling problems are considered in this 
study. Considering that P(10) = 12 <C 10! = 3628800, a complete enumeration 
approach is used for pairing and only the permutation of a sequence is focused. 

3 Discrete State Transition Algorithm 

In the case a solution to a specific optimization problem is described as a state, 
then the transformation to update the solution becomes a state transition. With- 
out loss of generality, the unified form of discrete state transition algorithm can 
be described as: 



( X k+1 = A fe (xfe)0Pfe(Mfe) 

\Vk+i = /(asfe+i) 



(7) 
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n 



Fig. 1. Growth trend relative to the sequence size 

where, x k S Z n stands for a current state, corresponding to a solution of a 
specific optimization problem; u k is a function of x k and historical states; ^lfc(-), 
Bk(-) are transformation operators, which are usually state transition matrixes; 
@ is an operation, which is admissible to operate on two states; and / is the 
cost function or evaluation function. 

The following three transformation operators are denned to permute current 
solution [TP] : 
(1) Swap Transformation 

x k+1 = A s k wap (m a )x k , (8) 

where, A s ^" ap 6 K™ x " is the swap transformation matrix, m a is the swap factor, 
a constant integer used to control the maximum number of positions to be ex- 
changed, while the positions are random. Fig. [2] shows an example of the swap 
transformation with m a = 2. 




123456 153426 



Fig. 2. Illustration of the swap transformation 
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(2) Shift Transformation 

x k+1 = A s k '" ft (m b )x k , (9) 

where, A^ 111 ^ G R nxn is the shift transformation matrix, nib is the shift factor, 
a constant integer used to control the maximum length of consecutive positions 
to be shifted. Note that both the selected position to be shifted after and po- 
sitions to be shifted are chosen randomly. Fig. [3] shows an example of the shift 
transformation with mj = 1. 




123456 124536 



Fig. 3. Illustration of the shift transformation 

(3) Symmetry Transformation 

x k+1 = A s k ym (m c )x kl (10) 

where, A s k ym £ W nxn is the symmetry transformation matrix, m c is the symme- 
try factor, a constant integer used to control the maximum length of subsequent 
positions as center. Note that both the component before the subsequent posi- 
tions and consecutive positions to be symmetrized are created randomly. Fig. 2] 
shows an example of the symmetry transformation with m c = 0. 




123456 154326 



Fig. 4. Illustration of symmetry transformation 



() 



A Multiobjective State Transition Algorithm for Single Machine Scheduling 



4 Pareto Archived Strategy based on DSTA 

In state transition algorithm, the times of transformation are called search en- 
forcement (SE); as a result, after each transformation operator, a candidate state 
set S is generated. 

4.1 Non-dominated sort 

We use a sorting approach similar to the fast-non-dominated-sort proposed in 
[j~2] . described as follows: 

1: for each s € S do 

2: n s <- 

3: for each t G S do 

4: if t -< s then 

5: n s ^— n s + 1 

6: end if 

7: end for 

8: end for 

where, n s is the domination count, representing the number of solutions domi- 
nating solution s. After the non-dominated sort, the state with the least count 
will be stored as incumbent best for the next transformation operator. 

4.2 Pareto Archived Strategy 

We adopt a simple Pareto archived strategy to select current best, as follows: 
1: for each A; e A do 
2: if best -< Ai then 
3: A^- A- Ai 

4: else if Ai -< best then 
5: A <- A 

6: else 

7: A^AlJfrest 
8: end if 
9: end for 

where, A is the archive keeping all non-donominated solutions. 

4.3 Pseudocodes of the proposed algorithm 

The core procedure of the proposed algorithm can be outlined in pseudocodes: 
1: repeat 

2: State <!— operator {best, SE, n) 

3: best <— update Jbest(best, SE, n, data) 

4: Paretoset 4— update -archive{Pareto set, best) 

5: until the maximum number of iterations is met 
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where, State is the state set; operator stands for the three transformation op- 
erators, which are carried out sequentially; update-best is corresponding to the 
non-dominated sort, and update jar chive corresponds to the Pareto archived 
strategy. The data is the known information (cost saving matrix, due dates and 
processing times) about a specific scheduling problem. 

5 Experimental Results 

In order to test the performance of the proposed multiobjective state transition 
algorithm, two typical examples are used for comparison. In the following exper- 
iments, SE = 20, m a = 2, mt = l,m c = are adopted for parameter settings. 
The maximum number of iterations for are 100 and 1000 respectively for the two 
examples. The known data for Example 1, 2 are given in Table Q] and Tabled 
Table [3] and Table HJ respectively, and the corresponding results can be found in 
Tableland Table |U It is worth to note that, the pairing methodology used with 
complete enumeration and Cuckoo Search (CS) is based on a greedy approach 
by first selecting the pair that produces the highest cost savings, and then re- 
peating the same procedure for the remaining set of pairs in the sequence [7]. 
We can find that for Example 1, STA obtained a solution which can dominate 
the optimal solutions by enumeration and CS. From both examples, it is easy to 
find that some additional optimal solutions are achieved by STA. 

Table 3. The cost savings matrix for 10 jobs having the same material 



Jobl Job2 Job3 Job4 Job5 Job6 Job7 Job8 Job9 JoblO 



Jobl 





2.73 


2.1 


2.16 


2.66 


3.6 


2.46 


2.7 


2.46 


2.8 


Job2 


2.73 





2 


1.6 


4.3 


3.69 


2.3 


3.5 


2.76 


3.6 


Job3 


2.1 


2 





1.4 


3.51 


3.33 


2.52 


3.68 


2.52 


2.46 


Job4 


2.16 


1.6 


1.4 





2.17 


2.32 


2.72 


3.04 


2.04 


2.97 


Job5 


2.66 


4.3 


3.51 


2.17 





3.6 


4.05 


4.41 


2.7 


2.64 


Job6 


3.6 


3.69 


3.33 


2.32 


3.6 





2.58 


4.7 


3.44 


2.94 


Job7 


2.46 


2.3 


2.52 


2.72 


4.05 


2.58 





2.6 


2.88 


2.82 


Job8 


2.7 


3.5 


3.68 


3.04 


4.41 


4.7 


2.6 





3.64 


3.57 


Job9 


2.46 


2.76 


2.52 


2.04 


2.7 


3.44 


2.88 


3.64 





3.76 


JoblO 


2.8 


3.6 


2.46 


2.97 


2.64 


2.94 


2.82 


3.57 


3.76 






6 Conclusion 

A multiobjective state transition algorithm is presented for a single machine job 
shop scheduling problem. In this paper, a complete enumeration approach is used 
for pairing the jobs in a fixed sequence. Compared with a greedy-based approach 
used with both the complete enumeration method and the CS, experimental 
results show the effectiveness of the proposed algorithm in obtaining the true 
set of all Pareto optimal solutions. 
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Table 4. Due dates and processing times for a set of 10 jobs 



Job 


Due Date(days) 


Processing Time(hours) 


Jobl 


11 


14:00 


Job2 


2 


18:00 


Job3 


13 


15:00 


Job4 


14 


8:20 


Job5 


11 


17:20 


Job6 


9 


16:00 


Job7 


4 


19:40 


Job8 


6 


23:20 


Job9 


10 


20:00 


JoblO 


10 


19:20 



Table 5. Comparison results for the set of jobs presented in Table Q] and Tabled 



Approach 


Optimal solutions 


T 


C 


Complete 


(2-5)-(l-4)-3 


13 


8.31 


Enumeration 


(5-2)-(l-4)-3 


13 


8.31 




(2-5)-(4-l)-3 


13 


8.31 




(2-4)-(5-l)-3 


15 


8.62 


CS0 


(2-5)-(l-4)-3 


13 


8.31 




(5-2)-(l-4)-3 


13 


8.31 




(2-5)-(4-l)-3 


13 


8.31 




(2-4)-(5-l)-3 


15 


8.62 


STA 


(5-2)-(l-4)-3 


13 


8.31 




(2-5)-(l-4)-3 


13 


8.31 




(2-5)-(4-l)-3 


13 


8.31 




(5-2)-(4-l)-3 


13 


8.31 




(2-4)-(5-l)-3 


15 


8.62 
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Table 6. Comparison results for the set of jobs presented in Table [3] and Table [4] 



Approach 


Optimal solutions 


T 


C 


Complete 


(5-7)-(2-6)-(l-3)-(4-10)-(8-9) 


39 


16.45 


Enumeration 


(5-7)- (2-6)- (1-3)- (4-8)- (10-9) 


40 


16.64 




(5-7)-(2-6)-(l-3)-(4-8)-(9-10) 


40 


16.64 




(5-7)-(2-6)-(l-4)-(3-8)-(10-9) 


41 


17.34 




(5-7)-(2-6)-(l-4)-(3-8)-(9-10) 


41 


17.34 




(5-2)-(7-4)-(6-l)-(3-8)-(10-9) 


43 


18.06 




(5-2)-(7-4)-(6-l)-(3-8)-(9-10) 


43 


18.06 




(2-5)-(7-4)-(6-l)-(3-8)-(10-9) 


43 


18.06 




(2-5)-(7-4)-(6-l)-(3-8)-(9-10) 


43 


18.06 


CS0 


2-(7-5)-(6-l)-3-(4-10)-(8-9) 


39 


14.26 




(5-7)-(2-6)-(l-3)-(4-8)-(9-10) 


40 


16.64 




(5-7)-(2-6)-(l-4)-(3-8)-(10-9) 


41 


17.34 




(2-5)-(7-4)-(6-l)-(3-8)-(10-9) 


43 


18.06 




(2-5)-(7-4)-(6-l)-(3-8)-(9-10) 


43 


18.06 




(5-2)-(7-4)-(6-l)-(3-8)-(10-9) 


43 


18.06 


STA 


(5-7)-(2-6)-(l-3)-(4-10)-(8-9) 


39 


16.45 




(5-7)-(2-6)-(l-3)-(4-10)-(9-8) 


39 


16.45 




(5-7)-(2-6)-(l-3)-(4-8)-(10-9) 


40 


16.64 




(5-7)-(2-6)-(l-3)-(4-8)-(9-lQ) 


40 


16.64 




(5-7)-(2-6)-(l-4)-(3-8)-(10-9) 


41 


17.34 




(5-7)-(2-6)-(l-4)-(3-8)-(9-10) 


41 


17.34 




(5-2)-(7-4)-(6-l)-(3-8)-(10-9) 


43 


18.06 




(5-2)-(7-4)-(6-l)-(3-8)-(9-10) 


43 


18.06 




(2-5)-(7-4)-(6-l)-(3-8)-(10-9) 


43 


18.06 




(2-5)-(7-4)-(6-l)-(3-8)-(9-10) 


43 


18.06 
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